TECHNICAL  NOTES 


One-Dimensional  Wave  Bottom  Boundary  Layer 
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Abstract:  Six  one-dimensional-vertical  wave  bottom  boundary  layer  models  are  analyzed  based  on  different  methods  for  estimating  the 
turbulent  eddy  viscosity:  Laminar,  linear,  parabolic,  k — one  equation  turbulence  closure,  k-e — two  equation  turbulence  closure,  and 
k-tx^ — two  equation  turbulence  closure.  Resultant  velocity  profiles,  bed  shear  stresses,  and  turbulent  kinetic  energy  are  compared  to 
laboratory  data  of  oscillatory  flow  over  smooth  and  rough  beds.  Bed  shear  stress  estimates  for  the  smooth  bed  case  were  most  closely 
predicted  by  the  k-m  model.  Normalized  errors  between  model  predictions  and  measurements  of  velocity  profiles  over  the  entire 
computational  domain  collected  at  15°  intervals  for  one-half  a  wave  cycle  show  that  overall  the  linear  model  was  most  accurate.  The  least 
accurate  were  the  laminar  and  k-e  models.  Normalized  errors  between  model  predictions  and  turbulence  kinetic  energy  profiles  showed 
that  the  k-(a  model  was  most  accurate.  Based  on  these  findings,  when  the  smallest  overall  velocity  profile  prediction  error  is  required,  the 
processing  requirements  and  error  analysis  suggest  that  the  linear  eddy  viscosity  model  is  adequate.  However,  if  accurate  estimates  of  bed 
shear  stress  and  TKE  are  required  then,  of  the  models  tested,  the  k-m  model  should  be  used. 
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Introduction 

Nearshore  surface  wave  motions  induce  flow  oscillations  near  the 
sea  bed  that  are  altered  by  frictional  resistance  in  what  is  termed 
the  wave  bottom  boundary  layer  (WBBL).  WBBLs  are  important 
for  a  variety  of  nearshore  and  coastal  engineering  problems,  in¬ 
cluding  sediment  transport  resulting  from  turbulence  and  turbu¬ 
lent  bed  shear  stresses  driven  by  these  oscillations.  Although 
much  is  known  about  steady  flow  boundary  layers,  study  of  os¬ 
cillatory  boundary  layers  still  gamers  much  research  mostly  due 
to  the  many  possibilities  for  parameterizing  the  turbulence. 

The  most  numerically  simple  methods  for  parameterizing  the 
turbulence  are  to  specify  the  shape  of  the  profile  a  priori  based  on 
the  friction  velocity  [e.g.,  Trowbridge  and  Madsen  (1984)  and 
Eredspe  and  Deigaard  (1992)]  whereas  the  computationally  ex¬ 


pensive  methods  require  direct  numerical  simulations  where  the 
model  grid  is  so  small  that  all  the  necessary  scales  of  motion  are 
calculated  rendering  a  parameterization  unnecessary  [e.g.,  Spalart 
(1988)].  A  middle  ground  approach  requires  descriptions  of  the 
turbulent  kinetic  energy  (TKE)  and  energy  dissipation  rate  in  the 
model  that  are  used  to  determine  the  turbulent  eddy  viscosity. 
These  turbulence  closure  schemes  (namely  k,  k-e,  and  k-ta) 
along  with  the  specified  eddy  viscosity  models  will  be  investi¬ 
gated  here  by  comparison  to  laboratory  data. 


Boundary  Layer  Theory  and  Numerical  Methods 

The  time-dependent  equation  of  turbulent  motion  for  an  incom¬ 
pressible  Newtonian  fluid  with  an  ensemble  mean  defect  oscilla¬ 
tory  velocity  component,  u^,  parallel  to  the  bed  is  given  by 
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where  v= kinematic  molecular  viscosity  of  the  fluid; 
r’,=turbulent  eddy  viscosity  in  the  usual  Reynold’s  averaged 
sense;  z=vertical  coordinate;  and  u^=u-U,  represents  the  differ¬ 
ence  between  the  depth-dependent  velocity  inside  the  boundary 
layer,  u  and  that  of  the  free  stream,  U. 

A  variety  of  methods  have  been  arrived  at  to  determine  or 
specify  v,.  One  option  is  to  assume  or  require  the  flow  be  laminar 
in  which  case  v,  is  zero  for  all  time.  Many  other  options  exist 
including  a  linear  form  using  a  time  varying  friction  velocity 
(Trowbridge  and  Madsen  1984),  a  parabolic  form  (Eredspe  and 
Deigaard  1992)  and  the  turbulent  closure  schemes  mentioned 
above.  We  adopt  the  one-equation  k  and  two-equation  k-e  model 


322  /  JOURNAL  OF  WATERWAY,  PORT,  COASTAL,  AND  OCEAN  ENGINEERING  ©  ASCE  /  NOVEMBER/DECEMBER  2004 


Report  Documentation  Page 

Form  Approved 

0MB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  0MB  control  number. 

1.  REPORT  DATE 

2QQ^  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2004  to  00-00-2004 

4.  TITLE  AND  SUBTITLE 

One-Dimensional  Wave  Bottom  Boundary  Layer  Model  Comparison: 
Specific  Eddy  Viscosity  and  Turbulence  Closure  Models 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Center  for  Applied  Coastal  Research  University  of 

Delaware, Newark, DE, 19716 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR’S  ACRONYM(S) 

11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  aS 

unclassified  unclassified  unclassified  Report  (SAR) 

4 

Standard  Form  298  (Rev.  8-98} 

Prescribed  by  ANSI  Std  Z39-18 


C- Parabolic  F-k-co 


Fig.  1.  Rough  bed  (Test  13),  experimental  data  (dots),  model  (solid 
curves):  (A)  Laminar,  (B)  linear,  (C)  parabolic  eddy  viscosity;  and 
(D)  k,  (E)  k-E,  and  (F)  k-w  turbulence  closure  schemes.  Each  zero 
velocity  tick  mark  is  offset  by  2  m  s“*. 


described  in  (Wilcox  2000)  and  the  k-u>  model  described  in  (Wil¬ 
cox  1988).  For  brevity  here,  complete  descriptions  of  the  models 
are  not  given  and  the  reader  is  referred  to  the  references  above  or 
Rodi  (1980),  Justesen  (1988),  or  Pope  (2000)  among  many  others 
for  further  description. 

The  equations  for  the  boundary  layer  flow  and  the  various 
eddy  viscosity  approaches  are  solved  numerically  via  an  implicit 
algorithm.  A  logarithmic  vertical  grid  is  used  with  asymmetrical 
second-order  spatial  differences.  Full  description  of  the  numerical 
technique,  boundary  conditions  and  forcing  mechanisms  are  pro¬ 
vided  in  Puleo  et  al.  (2003). 


Model  Results  and  Model-Data  Comparison 

The  WBBL  models  will  be  tested  against  Test  10  and  Test  13 
from  Jensen  et  al.  (1989).  Data  was  collected  at  every  15°  of 
wave  phase  for  eighty  wave  cycles  in  a  U-shaped  oscillatory 
water  tunnel  using  a  laser  Doppler  anemometer  (LDA).  Test  13 
was  forced  with  a  sinusoid  of  9.72  s,  a  maximum  velocity  of 
2  m  s“^  and  an  equivalent  Nikuradse  sand  grain  roughness  of 
Kf^=0.?i4  mm.  Test  10  was  forced  with  the  same  conditions  ex¬ 
cept  the  bed  was  made  smooth  by  the  installation  of  PVC  plates. 
Bed  shear  stress  was  also  measured  for  Test  10  using  a  hot  film 
probe.  Because  the  models  require  a  Zo  value,  a  smooth  bed  was 
approximated  by  using  1  X  10“'^  mm  (4  orders  of  magnitude 
smaller  than  the  rough  bed  case  and  one  order  of  magnitude 
smaller  than  clay  particle  size  on  the  Wentworth  scale).  Two  thou¬ 
sand  time  steps  per  wave  period  and  300  logarithmically  spaced 
grid  points  were  used  in  our  calculations. 

Velocity  Profiles 

Fig.  1  shows  the  velocity  comparison  for  every  15°  of  phase  for 
the  first  half  of  the  wave  for  the  simulations  (solid  lines)  and  the 
Jensen  et  al.  (1989)  Test  13  data  (dots).  Each  successive  profile  is 
offset  by  2  m  s“*  and  is  denoted  by  the  zeroes  in  the  plot.  Because 
of  the  turbulent  nature  of  the  flow,  the  laminar  (the  term  laminar 


Fig.  2.  Bed  shear  stress  for  (a)  Test  10;  smooth  bed  and  (b)  Test  13; 
rough  bed.  Laminar  (thin  dashed);  linear  (medium  dashed);  parabolic 
(bold  dashed);  k  (thin  solid);  k-s  (medium  solid);  and  k-w  (bold 
solid).  Test  10  data  (dots). 


here  refers  to  the  fact  that  only  molecular  kinematic  viscosity  is 
included  in  the  solution)  simulation  is  not  expected  to  match  the 
observed  profile  shape  and  is  given  for  illustrative  purposes. 

The  k-B  simulation  generally  has  the  correct  velocity  profile 
shape  but  fails  to  predict  the  magnitude  in  either  the  upper  portion 
of  the  domain  or  near  the  bed  where  it  over  or  under  predicts  the 
velocity  magnitude  at  nearly  all  wave  phases.  The  rest  of  the 
simulations  are  very  similar  and  qualitatively  match  the  data 
(quantitative  comparisons  are  given  in  the  following).  Smooth 
bed  simulation  (Test  10)  velocity  profile  comparisons  (not  shown) 
are  similar  to  those  of  Test  13  and  suggest  the  models  are  quali¬ 
tatively  similar  regardless  of  whether  or  not  they  are  applied  to 
smooth  or  rough  beds. 

Shear  Stress 

Bed  shear  stress  estimates  for  Tests  10  and  13,  shown  respec¬ 
tively,  in  Fig.  2(a  and  b),  indicate  variation  over  the  period  and 
between  the  different  model  formulations.  For  instance,  all  stress 
estimates  have  similar  shapes,  but  vary  in  magnitude  from  a 
maximum  of  L7kgm“*s“^  for  the  laminar  simulation  to 
17.4  kgm“'  s“^  for  the  k-s  model  simulations.  The  inflection  in 
the  k-B  model  bed  shear  stress  estimate  is  likely  attributed  to  the 
flow  reversal.  Bed  shear  stress  measurements  for  Test  10  are  most 
closely  predicted  by  the  k-w  model.  The  predictions  match  the 
data  for  most  of  the  wave  cycle  except  during  the  periods  of 
highest  flow  magnitude  and  capture  the  correct  timing  of  flow 
reversals.  The  k  simulation  also  captures  the  correct  timing  of 
flow  reversals,  but  over  predicts  the  bed  shear  stress  by  up  to 
factor  of  two.  The  linear  and  parabolic  eddy  viscosity  simulations 
over  predict  the  measured  bed  shear  stress  and  predict  the  minima 
in  bed  shear  stress  before  they  are  measured.  Finally,  the  k-B 
model  over  predicts  the  measurements  by  about  500%  and  is 
likely  the  cause  for  the  poor  predictive  capability  of  the  model  for 
this  case.  The  bed  shear  stress  estimates  for  Test  13  [Fig.  2(b)]  are 
similar  to  their  smooth  bed  counterparts  except  that  the  k-B 
model  estimate  overlays  the  L-co  model  estimate  further  suggest¬ 
ing  that  the  k-B  model  has  difficulty  on  the  smooth  bed  simula¬ 
tion. 
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k  (m^  s'^) 


Fig.  3.  Turbulent  kinetic  energy  for  Test  13  at  four  phases  for  k  (thin 
line),  k-E  (medium  line),  fc-co  (thick  line),  and  data  (dots) 


Turbulent  Kinetic  Energy 

Fig.  3  compares  TKE  profiles  at  four  phases  of  the  wave  cycle  for 
the  k,  k-B,  and  k-ia  models  to  the  data  for  Test  13.  TKE  was 
estimated  from  laboratory  measurements  as  0.65  times  the  sum  of 
the  mean  square  velocity  fluctuations  in  the  flow-parallel  and 
flow-normal  directions  (Justesen  1991).  It  is  clear  from  the  simu¬ 
lations  that  the  k-B  model  has  difficulty  reproducing  the  mea¬ 
sured  TKE  near  the  bed  regardless  of  the  phase  of  the  free  surface 
flow.  In  other  words,  our  implementation  of  the  k-s  model  shows 
difficulty  in  predicting  the  TKE  even  away  from  strong  adverse 
pressure  gradients.  In  contrast,  the  k  and  k-d)  models  more 
closely  match  the  measurements,  even  though  variations  are  seen 
both  near  the  bed  and  in  the  upper  portion  of  the  boundary  layer. 
Similar  results  were  found  for  Test  10,  but  for  brevity  are  not 
shown  here.  Overall  prediction  errors  for  both  tests  are  described 
below. 


Eddy  Viscosity  Models  versus  Turbulence  Closure 
Schemes 


A  quantitative  measure  of  the  predictive  capability  of  the  models 
can  be  determined  from  the  normalized  phase  dependent  error 
(the  sum  over  z  of  the  estimated  error  variance  normalized  by  the 
variance  in  the  data) 


p  ^  V '  (Xdata  Xmodel) 


(2) 


where  crj^jj|(tp)= velocity  variance  (over  the  vertical  direction)  in 
the  measurement  data  at  each  phase  and  x  model  prediction;  in 
this  study  either  the  velocity  or  TKE.  The  smooth  and  rough  bed 
comparisons  for  velocity  show  that  the  linear  model  was  most 
accurate  (although  the  parabolic,  k  and  k-u>  models  had  similar 
errors)  whereas  the  laminar  model  was  least  accurate  (Fig.  4; 
laminar  errors  extend  beyond  the  axis  range  shown).  The  k-B 
model  had  normalized  errors  that  were  up  to  6  times  as  large  as 
the  k  and  fc-w  models.  The  errors  for  velocity  become  largest 
between  0  and  20°  and  between  160  and  180°  when  the  flow  is 
decelerating  and  passing  through  flow  reversal.  Away  from  flow 
reversals  (adverse  pressure  gradients),  all  of  the  models  except 
the  laminar  and  k-B  perform  nearly  equally  for  both  the  smooth 
bed  [Fig.  4(a)]  and  rough  bed  [Fig.  4(b)]  tests.  In  general  the  best 


Fig.  4.  Normalized  phase  dependent  velocity  errors:  (A)  smooth  bed 
Test  10;  and  (B)  rough  bed  Test  13.  Normalized  turbulent  kinetic 
energy  (TKE)  errors  for  the  turbulent  closure  schemes:  (C)  smooth 
bed  Test  10  and  (D)  rough  bed  Test  13.  Laminar  (gray  asterisk  dash); 
linear  (asterisk);  parabolic  (square);  k  (diamond);  k-E  (circle);  and 
k-m  (cross). 


performers  for  velocity  were  the  linear  and  parabolic  eddy  viscos¬ 
ity  simulations  and  the  X-u  model. 

Figs.  4(c  and  d)  show  the  normalized  TKE  errors  but  for  the  k, 
k-B,  and  k-ia  models  only.  Again,  it  is  clear  that  the  k  and  k-ia 
models  outperform  the  k-B  model  especially  in  the  smooth  bed 
simulation  [Fig.  4(c)]  case  where  the  errors  for  the  k-B  model 
extend  well  outside  the  axis  range. 

The  overall  predictive  capability  of  the  models  can  be  deter¬ 
mined  from  the  normalized  total  error  (normalized  error  variance 
summed  over  space  and  summed  over  each  phase). 

This  calculation  was  performed  for  the  velocity  only  for  the 
linear  and  k-ta  model  for  time  steps  ranging  from  100  to  2,000 
per  wave  cycle  to  address  the  predictive  error  as  a  function  of 
computation  time.  It  was  found  that  the  velocity  errors  were  not 
statistically  different  for  any  of  the  various  time  steps.  However, 
the  fastest  k-ta  simulation  was  1.5  times  slower  than  the  slowest 
linear  eddy  viscosity  simulation  and  about  40  times  slower  than 
the  fastest  linear  eddy  viscosity  simulation.  The  same  calculations 
were  not  carried  out  for  TKE  due  to  difficulty  in  accurately  esti¬ 
mating  TKE  from  the  linear  eddy  viscosity  model. 


Conclusions 

Six  one-dimensional  wave  bottom  boundary  layer  models  have 
been  compared  to  laboratory  data  on  smooth  and  rough  beds.  In 
general,  the  eddy  viscosity  models,  where  the  shape  of  the  verti¬ 
cal  profile  of  eddy  viscosity  is  specified  perform  just  as  well  or 
better  than  one  and  two  equation  turbulence  closure  schemes  in 
predicting  the  overall  velocity  profile.  In  addition,  it  was  found 
that  the  model  using  a  linear  profile  of  eddy  viscosity  yields 
slightly  better  results  than  the  more  numerically  intensive  X-co 
model  on  a  rough  bed  at  about  40  times  savings  in  computer 
processing  requirements.  However,  the  k  and  k-ia  models  had  a 
relatively  high  predictive  capability  in  terms  of  TKE  and  the  k-ia 
model  was  the  most  accurate  predictor  of  bed  stress.  It  must  be 
kept  in  mind  that  the  findings  in  this  note  are  for  high  Reynolds 
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number  (>10^)  flows  and  have  not  been  verified  or  tested  for 
small  Reynolds  number  flows  in  this  study.  In  addition,  these 
findings  are  not  easily  extrapolated  to  multidimensional  (2D  or 
3D)  models,  but  suggest  that  for  one-dimensional  flat  bed 
(smooth  or  rough)  simulations  no  increase  in  predictive  capability 
is  obtained  in  using  the  k,  k-e,  or  k-m  turbulence  closure 
schemes  over  the  more  simplified  linear  or  parabolic  eddy  viscos¬ 
ity  models  with  respect  to  total  error  variance  for  velocity  pro¬ 
files.  However,  if  one  is  interested  in  accurately  predicting  TKE 
or  more  importantly  for  sediment  transport  studies,  the  bed  shear 
stress,  then,  of  the  models  tested,  the  fc-co  model  is  recom¬ 
mended. 
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Notation 

The  following  symbols  are  used  in  this  technical  note: 

Ep  =  normalized  phase  dependent  error  between  model  and 
data; 

e  =  turbulent  dissipation  rate; 

Kff  =  Nikuradse  equivalent  sand  roughness; 
k  =  turbulent  kinetic  energy; 

=  number  of  grid  points; 
t  =  time; 


U  =  free  stream  fluid  velocity; 
u  =  bed  parallel  fluid  velocity; 

=  deflect  velocity  (m-U); 

Zq  =  zero  level  for  velocity  (kN/30); 
z  =  distance  above  horizontal  bed; 

V  =  kinematic  viscosity; 

V,  =  turbulent  eddy  viscosity; 
cr^  =  velocity  variance; 
cp  =  oscillatory  flow  phase; 

X  =  model  prediction  for  velocity  or  TKE;  and 
CO  =  specific  dissipation  rate  {elk). 
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